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Abstract 



Gravitational wave experiments will play a key role in the investigation 
of the frontiers of cosmology and the structure of fundamental fields at high 



, energies, by either setting stringent upper limits on, or by detecting the pri- 

mordial gravitational wave background produced in the early-Universe. Here 
we discuss the impact of space-borne laser interferometric detectors operating 
C*~) ■ in the low-frequency window ~ 10~ 6 Hz - 1 Hz; the aim of our analysis is 

to investigate whether a primordial background characterized by a fractional 
energy density /i 2 00 ^ ~ 1CP 16 - 1CP 15 , which is consistent with the prediction 
*q ' of "slow-roll" inflationary models, might be detectable by the Laser Interfer- 

ometer Space Antenna (LISA) or follow-up missions. 

In searching for stochastic backgrounds, the presently planned LISA mis- 
sion suffers from the lack of two detectors with uncorrelated noise. We analyze 
the sensitivity improvements that could be achieved by cross-correlating the 
data streams from a pair of detectors of the LISA class; we show that this 
configuration is extremely compelling, leading to the detection of a stochas- 
tic background as weak as /i 2 00 $7 ~ 5 x 10~ 14 . However, such instrumental 
sensitivity can not be fully exploited to measure the primordial component 
of the background, due to the overwhelming power of the signal produced 
by large populations of short-period solar-mass binary systems of compact 
objects. We estimate that the primordial background can be observed only if 
its fractional energy density /i 2 00 VL is greater than ~ 5 x 10~ 13 . 

The key conclusion of our analysis is that the stochastic radiation from 
unresolved binary systems sets a fundamental limit on the sensitivity that 
can be achieved in searching for the primordial background in frequencies 
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between ~ 10 Hz and 0.1 Hz, regardless of the instrumental noise level and 
the integration time. Indeed, the mHz frequency band, where LISA achieves 
optimal sensitivity, is not suitable to probe slow-roll inflationary models. We 
briefly discuss possible follow-up missions aimed at the frequency region ~ 
0.1 Hz - 1 Hz, which is likely to be free from stochastic backgrounds of 
astrophysical origin: no fundamental limits seem to prevent us from reaching 
/if 00 $7 ~ 10~ 16 , although the technological challenges are considerable and 
deserve careful study. 



PACS numbers: 04.30.-y, 04.80.Nm,98.80.Es 
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I. INTRODUCTION 



The Universe became "thin" to gravitational waves (GWs) at the Plank epoch, corre- 
sponding to the cosmic time ~ 10 -43 sec; the gravitons decoupled from the surrounding 
plasma at a temperature of the order of the Planck mass ~ 10 19 GeV, and gravitational 
radiation produced at that epoch or later - including the electro-weak and the Grand Uni- 
fication (GUT) scale - has travelled undisturbed to us, carrying full information about the 
state of the Universe, and the physical processes from which it took origin. Indeed GW ex- 
periments will open radically new frontiers for cosmology and high energy physics (see 
and reference therein for an extensive discussion). 

In the time frame ~ 2002-2010 a large portion of the GW spectrum will progressively 
become accessible, mainly through large-scale laser interferometers. On the ground, the 
world-wide network of km-size interferometers - LIGO, GEO600, VIRGO and TAMA - sen- 
sitive in the frequency band ~ 10 Hz — 1kHz, will start carrying out "science runs" at the 
beginning of 2002, with the realistic goal of directly detecting GWs. Several instrumental 
upgrades, starting around 2005, will drive the sensitivity of the instruments to a GW stochas- 
tic background from hf 00 fl ~ 10 -6 (for the so-called initial generation) to h\ QQ Q ~ 10~ 10 (for 
the so-called advanced configuration). In space, a collaboration between ESA and NASA 
is carrying out the project called LISA (Laser Interferometer Space Antenna). This is a 
space-borne laser interferometer with arms of length 5 x 10 6 km, planned to fly by 2010 [[3]. 
This instrument guarantees the detection of GWs at low frequencies (~ 10~ 5 Hz — 10~ 2 Hz). 

The purpose of this paper is to show the central role of the experiments in the low- 
frequency window ~ 10~ 6 Hz - 1 Hz, with emphasis on instruments of the LISA class, in the 
search for the primordial GW background. Our aim is to identify the fundamental issues 
regarding the achievement of a sensitivity in the range hf 00 Q ~ 10~ 16 - 10 -15 , which is set 
by the theoretical prediction of "slow-roll" inflationary models. 



A. The stochastic background spectrum 

A stochastic GW background is a random process that can be described only in terms 
of its statistical properties. Without loss of generality, for the issues discussed in this paper, 
we assume it to be isotropic, stationary, Gaussian and unpolarized. The energy and spectral 
content of a stochastic background are described by the dimensionless function 

n(/)s Ife(fl (L1) 
Pc a In / 

p gw is the energy density carried by the background radiation, and 

Pc = T^TT- ~ L6 x 10_8 ^100 erg/cm 3 , 

^1.2xlO- m h 2 100 sec- 2 (1.2) 

is the critical energy density required today to close the Universe. The value of the Hubble 
constant (today) is 
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H = 100 h wo km sec" 1 MpcT 1 , 

~ 3.2 x HT 18 h 100 sec" 1 , (1.3) 

where h wo is known from observations to be in the range 0.4 < h 100 < 0.85. fi(/) is therefore 
the ratio of the GW energy density to the critical energy density per unit logarithmic 
frequency interval; one usually refers to hf 00 Q(f), which is independent of the unknown 
value of the Hubble constant. 

It is useful to introduce the characteristic amplitude h c (f) of the GW background: it is 
the dimensionless characteristic value of the total GW background-induced fluctuation h(t) 
at the output of an interferometer per unit logarithmic frequency interval: 



OO 



(h 2 (t))=2 d(\n f)hi(f) ; (1.4) 

JO 

here ( ) denotes the expectation value. The spectral density S(f) of the background is related 
to h c (f) by 

h 2 c U) = 2fS(f) , (1.5) 
and fi(/), h c (f), and S(f) satisfy the relation [l|] 

Q{f) = ^m fhlU) = M fsU) - (L6) 

The characteristic amplitude over a frequency band Af is therefore: 
h c (f,Af) = h c (f) (^pj ' 



7.1 x 10 -22 



10- 8 



1/2 / , v-3/2 / vl/2 



lmHz \3.2xlO- 8 Hz 



;i.7) 



where A^f ~ 3.2 x 10 8 (lyr/T) Hz is the width of the frequency bin for an observation 
time T. For comparison, the relevant characteristic amplitude of the LISA noise is ~ 10~ 24 . 



B. Sources of stochastic backgrounds 

The stochastic GW background can be divided into two broad classes, based on its 
origin: (i) the primordial GW background (PGB), produced by physical processes in the 
early Universe, and (ii) the astrophysically- generated GW background (GGB), generated 
by the incoherent superposition of gravitational radiation produced, at much later cosmic 
times, by a large number of astrophysical sources that can not be resolved individually. The 
emphasis of this paper is on the detectability of the PGB. 

In this paper we will use the following conventions: Q p (f) and £l g (f) identify the frac- 
tional energy density in GWs, Eq. (|1 . 1| ) , carried by the primordial and the generated com- 
ponent of the GW background, respectively. If no index is used, we refer to a general GW 
stochastic signal, with no assumption about its production mechanism. 

At present, there are three observational constraints on the PGB contribution to f2(/): 
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1. 



The high degree of isotropy of the cosmic microwave background radiation sets a limit 
at ultra- low frequencies @]: 




h 2 10Q n p (f) < 7 x 10- 11 (-L) , 3 x 1(T 18 h m Hz < / < 10- 16 h 100 Hz ; (1.8) 



2. The very accurate timing of milli-second radio-pulsars constrains Q p (f) in a frequency 
range of the order of the inverse of the observation time, typically of order of a few 
years [|]: 

hlMf) < 1(T 8 /~l(T 8 Hz; (1.9) 

3. The standard model of big-bang nucleosynthesis constrains the total energy content 
in GWs over a wide frequency range ||: 

/ ^ 00 fi p (/Mln/)<6xl(r 6 . (1.10) 

To foresee what physical processes could have produced a detectable GW background is 
an almost impossible challenge; nonetheless, it is enlighting to discuss some general principles 
and possible generation mechanisms to show the typical sensitivity that experiments should 
achieve in order to test different models. 

The main mechanisms that produce a PGB can be divided into two broad categories (for 
a recent detailed review see [0): (i) Parametric amplifications of metric tensor perturbations 
that occurred during an inflationary epoch, and (ii) Some causal processes - such as phase 
transitions - that took place in the early Universe. 

Stochastic backgrounds produced by the parametric amplification of metric tensor per- 
turbations that occurs during an inflationary epoch J7| extend over a huge range of fre- 
quencies, from ~ 3 x 10~ 18 Hz up to a cutoff frequency in the GHz range. In the win- 
dow ~ 10~ 16 Hz — 1 GHz, slow-roll inflationary models predict a quasi scale-invariant spec- 
trum whose typical magnitude - in order to satisfy the COBE bound - cannot exceed 
^ioo ~ 10~ 14 in the LISA frequency band, as well as in the Earth-based detectors ob- 
servational window ||; a more refined analysis [[J yields a more conservative upper limit: 
hf 0Q Q p ~ 10~ 16 — 10~ 15 . Superstring-inspired cosmological models |lO|-|T2l| predict a spectrum 



that, for suitable choices of the free parameters of the model, could reach hf Q0 Q p ~ 10~' 
at the frequencies accessible either to Earth-based or to space-borne experiments, while 
satisfying the existing observational bound p3HTT| . 

Stochastic backgrounds can also be produced by some classical causal processes that took 
place in the early Universe; for this class of signals, the characteristic frequency is related 
both to the time of emission and the corresponding temperature T. 

Non-equilibrium processes that occur at the reheating that takes place after inflation 
could provide a stochastic background with cutoff frequency in the range ~ 0.1 mHz — 1 kHz, 
corresponding to reheating temperatures between ~ 1 TeV and ~ 10 9 GeV. As an example, 
in hybrid and extended inflationary models, the exit towards a radiation-dominated era is 
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characterized by a first-order phase transition, which, if strongly of the first order, generates 
a stochastic background with h\ m Q p ~ at frequencies that can vary from the LISA 



observational window up to the sensitivity band of Earth-based interferometers JT! 

Phase transitions that inevitably occur at T ~ 10 2 MeV (the QCD phase transition) 
and T ~ 10 2 GeV (the electroweak phase transition) produce GWs. In particular, if the 
electroweak phase transition is strongly of the first order, the spectrum is approximately 
hf Q0 Q p ~ 10~ n — 10~ 9 at / ~ lmHz [1BJ; the requirement of a strong first order phase 



"100 J "V 

transition, which is necessary in order to have baryogenesis at the electroweak scale (see 



and references therein for a recent review), is directly related, in a minimal supersymmetric 



extension of the standard model, to the mass of the super-partner of the top quark [21-23 



Cosmic strings, which are topological defects formed during phase transitions, produce 
GWs whose typical frequency ranges from / ~ 10~ 8 Hz up to / ~ 10 10 Hz with /i 2 00 Q p ~ 
10 -9 — 10~ 8 , see |24| and references therein for a review. 

Global phase transitions associated with some scalar field which acquires a non-zero 
vacuum expectation value (VEV) below a critical temperature would produce, via a quite 
general relaxation process, GWs whose energy content is very significant, /i 2 00 Q p ~ 10~ 6 , 



for VEV's near the Planck/string scale p5| , |26 



Recently there has been a great amount of theoretical activity investigating higher di- 
mensional " brane- world" scenarios, where gravity begins to probe the extra dimensions at 
energies as low as 10 3 GeV; and an estimate of possible GW backgrounds in such models 
was presented recently in [|27| . 



These examples clearly show that the investigation of the primordial GW stochastic 
background in the low-frequency regime would provide us key information about the physics 
beyond the standard model and/or could allow us to discriminate between different infla- 
tionary cosmological models. 



C. Detecting a stochastic background 

A stochastic background is a random process which is intrinsically indistinguishable from 
the detector noise. In order to detect such a signal, the optimal signal processing strategy 
calls for correlations between two (or more pairs of) instruments, possibly widely separated 
in order to minimize the effects of common noise sources. The relevant data analysis issues 



have been thoroughly addressed in [^g,^J; here we simply review the main concepts, and 
refer to [[2^,[29[], and references therein, for more details. 

The statistical analysis is based on the following assumptions: the signal and the detector 
noise are uncorrelated; the noise in each detector is stationary and Gaussian, and possible 
noise correlations between two detectors are negligible. 

We define the output (signal + noise) of the two instruments as oi(t) and 02(f); the 
cross-correlation signal C that one constructs is therefore of the form: 

r T/2 rT/2 



C 



/ dt dt'o 1 {t)o 2 (t')Q(t-t'), (1.11) 

J-T/2 J-T/2 

where Q(t — t') is a suitable filter function. In the general case, the filter function depends 
on t and t' independently, that is Q = Q(t, t')\ here we have used the property of the signal 
of being stationary, and therefore Q(t,t') = Q(t — £'). The SNR is defined as: 
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SNR 



a 



where /x and a are the mean value and the variance of the observable C: 

' ml 



ti = {C) = T 



20vr 



1.121 



;i.i3) 



a* = (C*) - (Of = % (Q,Q) . 



;i.l4) 



;i.i5) 



where Q(f) is the Fourier transform of Q(t — t'). The functions R(f) and A(f) are defined 
as follows: 



Eqs. ( |1.13| ) and (|1.14| ) are written in terms of the usual inner product [p9| 

(0,6)= / dfa*(f)b(f)R(f), 



1 + 



IOtt 2 



/ 3 



+ 



3iJ 2 Vfi 2 (/) [l + 7 (/) 2 ] 



^(/)^ 2) (/) J V 107r V f*sP(f)Sk 2) {f) J ' 



1.16) 



tC/W) 
/ 3 W) ' 



;i.i7) 



In Eq. (P^D , S^(f),k = 1,2 is the one-sided noise power spectral density of the fc— th 
detector, and 7(7) is the so-called overlap reduction function, which depends entirely on the 
relative orientation and location of the two detectors; it accounts for SNR losses that occur 
when the instruments are not optimally located and oriented, cf. Eq. ( [L.20| ) and Sec. [TJ. 
Using Eqs. (|1.13|) and ( |1.14| ), one can cast Eq. ( |1.12|) in the form: 



SNR 2 



3HZ\ 2 (Q,£* 



IOtt 2 



Q,Q 



1.18) 



The optimal choice of the filter Q, is thus based on the maximizing the SNR, Eq. ( |1.18|) 
and is given by: 



Q(f) = (const.) x A(f) 



1.19) 



where the overall normalization factor is arbitrary. Note that Eqs. ( |1.11| )-( |I.19| ) are valid for 
a background of arbitrary energy density £l(f). In the case of a signal much weaker than 
the noise, i?gfi(/)// 3 <C S^(f), one can Taylor expand Eqs. ( |1.16| ) and ( |1.17| ), retaining 
only the leading order term. As a consequence, Eq. (|1.18| ) reduces to: 
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SNR 



50tt 2 



df 



o ' f*s£\f)Sk'>{f) 



7(2)/ 



1/2 



(signal <C noise) 



1.20) 



It is convenient to introduce the noise characteristic amplitude h Tms , equivalent to h c , as 
follows: 



dfS n (f) 
d(lnf)h 2 Ims (f). 



It is enlightening to write Eq. ( |1.20| ), using Eqs. (|1.4j) and ( 1.21 ), in the form: 



SNR~ 7 (/ C ) (A/T) 



1/2 



fec(/c) 
^rms(Jc) 



;i.2i) 



1.22) 



here we have assumed that the frequency band A/, which contains most of the SNR, is 
centered on the characteristic frequency f c , and is sufficiently small; the noise spectral 
density of the two instruments, that for simplicity we assume identical, and the overlap 
reduction function can be therefore treated as roughly constant. If only one instrument 
is in operation, one could in principle detect a stochastic background with SNR £ 1 when 
with two instruments one can detect the signal when h c £ h rms /[y(f c ) (A/T) 1 / 4 ]. 
Cross-correlation experiments are therefore highly desirable for both detection confidence 
and sensitivity. In fact, one can isolate the stochastic GW signal from all the spurious 
contributions which are uncorrelated between the two instruments. Common noise sources, 
which correlate on the same light-travel time scale, might, however, be present, degrading 
the overall sensitivity. Moreover, GW signals are expected to be very weak, well buried into 
the noise; using cross-correlations, through optimal filtering one increases the sensitivity by 
a factor ~ 10 (A//1 mHz) 1//4 (T/10 7 sec) 1//4 , with respect to the single detector case. 



D. Summary of the results 

The goal of this paper is to explore the capability of space-borne laser interferometers, 
such as LISA and its successors, in searching for the primordial GW stochastic backgrounds. 
The analysis of the LISA technology leads us to the following conclusions: 

• A PGB of fractional energy density h 2 00 Q p £ 10~ 10 definitely shows up as an excess 
power component in the data of a single LISA interferometer over a large frequency 
window; it might be detectable by calibrating the noise-only response of the instru- 
ment, but the issue of decisively assigning this contribution to a real primordial signal 
remains open. 

• Cross-correlations between the data streams of two identical LISAs, characterized by 
the presently estimated instrumental noise, allow us to reach a minimum value of the 
fractional energy density carried by GWs in the range 5 x 10~ 14 £ /i^ 00 fi^ mm ^ £ 10~ 12 
for an integration time T = 10 7 sec, depending on the location and orientation of the 
two detectors. 
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• Such remarkable sensitivity, however, does not apply to primordial GW backgrounds; 
in fact, the copious number of short period solar- mass binary systems in the Universe 
produce a GGB that overwhelms the PGB in the key mHz region; we estimate that the 
minimum detectable value of the primordial GW background is h\ 00 f2^ mm ) £ 5 x 1CT 13 . 

• The cross-correlation of the data streams from two LISA detectors provides, therefore, 
a powerful tool to extract information about populations of binary systems of short- 
period solar-mass compact objects in the Universe: the GGB is detectable at signal- 
to-noise ratio ~ 100. 

We would like to emphasize that the third generation Earth-based laser interferometers, 
instruments such as advanced LIGO (LIGO III) and EUGO, a new European detector 
currently under study, will be able to achieve a sensitivity /i 2 00 fi( mm ) ~ 10~ 10 ; space-based 
interferometers will therefore play a primary role in GW cosmology. 

Although the former results are very encouraging, our analysis leads to the rather 
obvious, but somewhat disappointing outcome that experiments in the frequency band 
~ 10~ 6 Hz — 0.1 Hz will be limited to a sensitivity of the order h\ QQ Q p ~ 10~ 13 — 10~ 12 ; 
this limit can not be improved by reducing the instrumental noise and/or increasing the 
integration time; in fact, GGB's produce a residual correlation in the filter output designed 
to detect the PGB, that can not be eliminated: the mHz frequency window is not suitable to 
search for a primordial gravitational wave background characterized by h\ 00 Q p < 5 x 10~ 13 . 

It is therefore worth asking whether future technological developments, and more am- 
bitious missions will enable the detection of a very weak PGB. Our present understanding 
suggests that cross-correlation experiments carried out with a pair of space-based interfer- 
ometers with optimal sensitivity in the band ~ 0.1 Hz — 1 Hz could be able to meet the target 
ft-ioo ~ 10~ 16 . In this frequency band, in fact, astrophysically generated backgrounds are 
not present, and an experiment is limited in principle only by the instrumental noise. The 
design of a detector aimed at the window ~ 0.1 Hz — 1 Hz imposes stringent requirements on 
the power and frequency of the laser, as well as on the dimensions of the "optics" and on sev- 
eral other components of the instrument. In order to test slow-roll inflation an effective - i.e. 
after residual radiation from individual binary systems has been removed - rms noise level 
~ 10~ 24 , and rather long integration times (~ 3 yrs) are required. Such technological and 
data analysis issues have been little investigated so far, and deserve careful consideration. 

The paper is organized as follows. In Sec. |TI| we derive a closed form expression of 
the overlap reduction function for a pair of space-based interferometers. In Sec. |Hj we 
review our present astrophysical understanding of the GW background generated by the 
incoherent superposition of radiation emitted by galactic and extra-galactic short period 
solar-mass binary systems. In Sec. [IV] we present the key results of the paper: we estimate 
the sensitivity that can be achieved by the present LISA technology, and possible follow-up 
missions, by cross-correlating the outputs of two identical detectors with uncorrelated noise. 
Sec. [V] contains our conclusions and pointers to future work. 
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II. THE OVERLAP REDUCTION FUNCTION FOR LISA DETECTORS 



In this section we compute a closed form expression for the overlap reduction function 
7(/), see Eqs. ( |L16| )-( |ITl~9| ), for space-borne interferometers, and in particular for instru- 



ments characterized by a LISA-like orbital configuration. Indeed, our analysis provides 
explicit formulae that can be directly applied, with little changes, to any orbital configura- 
tion. 

The overlap reduction function is a dimensionless function of the frequency /, which 
measures the degradation of the SNR when the detectors are not optimally oriented and 
located. At any given frequency, 7(/) depends entirely on the relative separation and ori- 
entation of the instruments, and its behavior has a very intuitive physical explanation. It 
is maximum when the detectors are co-located and co-aligned. As the detectors (paral- 
lel to each others) are shifted apart, at a distance D, the signal drives oscillations in the 
two instruments that are progressively out of phase; there is an effective correlation only 
if the separation is smaller than approximately half of the characteristic wavelength of the 
gravitational radiation; equivalently, the frequency band over which one accumulates SNR 
is 

/£1 (y^j) lmHz > (2-1) 

where one Astronomical Unit (AU) corresponds to 1.4959787 x 10 13 cm. For two coincident 
detectors, ■y(f) decreases, at any frequency, as one instrument is rotated with respect to the 
other, because the detectors are excited in different ways by the different polarizations; for a 
rotation of 7r/4, the overlap reduction function is identically zero. In the general case, 7(/) 
shows a complex behavior which is the superposition of the two effects that we have just 
described. 



A. The LISA mission 

LISA is an all-sky monitor with a quadrupolar antenna pattern. Its orbital configuration 
was conceived in order to keep the geometry of the interferometer as stable as possible during 
the mission, as well as to provide an optimal coverage of the sky: a constellation of three 
drag-free spacecraft (containing the free-falling test masses) is placed at the vertices of an 
ideal equilateral triangle with sides ~5x 10 6 km; it forms a three-arm interferometer, with 
a 60° angle between two adjacent laser beams. The LISA orbital motion is such that the 
barycentre of the instrument is inserted in a heliocentric (essentially circular) orbit, following 
by 20° the Earth; the detector plane is tilted by 60° with respect to the Ecliptic and the 
instrument counter-rotates around the normal to the detector plane with the same period 
of 1 yr. 

We introduce a Cartesian reference frame (x, y, z) tied to the Ecliptic, with z perpen- 
dicular to the Ecliptic plane, and x and y in the plane itself, and oriented in such a way to 
define a left-hand tern. In this frame, LISA's centre-of-mass is described by the polar angles 

® = \ , =$ + n e t , U EE^-j; (2.2) 
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$0 sets the position of the detector barycentre at some arbitrary reference time. The time 
evolution of the unit vectors lj (j = 1,2,3) along each arm are described by the following 
expression [50] : 



+ 



sinaj(i) cos$(t) — cos a j(t) sin$(i) 



sin a j(t) sin$(i) +cosot,-(t) cos$(£) 



y + 



V3 



sin a_j (t) 



where o.j{t) increases linearly with time, according to 



aj (t) = n®t - {j - 1)tt/3 + a , 



(2.3) 



(2.4) 



and «o is just a constant specifying the orientation of the arms at the arbitrary reference 
time t = 0. 

In the next section we will derive the expression of the overlap reduction function for 
two detectors characterized by a LISA-like motion. The time dependence of lj and the 



centre-of-mass are described by Eqs. Q2.2| ), (|2.3| ) and (|2.4j ) for both instruments, just with 
different initial conditions and <3>o- We will use the notation aoi and $oi, and CK02 and $02 
to indicate the corresponding values of the detector "1" and "2", respectively. For future 
convenience, we also introduce the notation: 



A$ = $ 



02 



01 



A« = a 2 - «01 • 



(2.5) 



The main noise sources that affect the mission have been addressed by the LISA Science 
Team and yield the following expression for the expected noise spectral density 



S n (f) ~ 1(T 48 1.22 x 1(T 3 f- 4 + 16.74 + 9.7 x 10 4 / 



4 r2 



Hz" 1 



(2.6) 



the rms noise amplitude that is inferred from S n (f) is shown in Fig. III. 



B. The overlap reduction function 



To derive a closed form expression of the overlap reduction function 7(/), we follow the 
formalism developed by Allen and Romano p9"| , which in turn was based on the analysis 
done by Flanagan [pjf| . j(f) is formally given by 



1(f) = Pi{x) d { ^d {2)ab + (*{x) d^S b d^ ac S c + p 3 (x) d%d$ S «S b S c S d . 
In the previous expression, are the detector response tensors, defined as: 



d 



(fc) 
06 



1 r 



(2.7) 



(2.8) 



where k = 1,2 labels the instrument, and rh^ and are the unit vectors along the arms 
of the interferometers; our (arbitrary) choice corresponds to li and I2, respectively, given 
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FIG. 1. The characteristic amplitude of stochastic gravitational wave backgrounds of primor- 
dial and astrophysical origin and the LISA rms noise amplitude. The plot shows the LISA rms 
noise amplitude in one year of observation (solid line), as a function of frequency, and compares 
it to the characteristic amplitude of several stochastic signals: a backgrounds with flat spectrum 
hi 00 O = 10~ 10 , 10~ 13 , 10~ 16 (dotted-lines, from top to bottom, respectively), and the astrophysi- 
cally generated galactic background (dashed- line) , according to the estimate given by Bender and 
collaborators, cfr. Eq. ([O]). 



by Eqs. Q2.3| ) and (|2.4|), with the appropriate initial conditions. The unit vector along the 
direction that connects the centers of mass of the two detectors is 

~ cos $9 — cos $1 sin $9 — sin $i 

" " (2.9) 



2(l-cos($ 2 -$ 1 )) J2(l - cos($ 2 - $0) 
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x is a dimensionless parameter defined as 

x = —^, 2.10 

c 

where 

D = R^2 (l-cosA$ ), (2.11) 
and the functions Pj{x) are given by 

10 5 

pi(x) = 5j (x) ji(x) + -^j 2 (x) , 

X X A 
/ N -l r\ • / \ 40 . . N 50 . , N 

P2O) = -lOjoW "I JiW 2 > 

5 25 175 

P 3 (x) = - Jo(^) - — ji(x) + ^J2(x) , (2.12) 



where 



sin x 

Jo Or) 



Ji(^) 



x 

sin x cos x 



x 2 x 



smx cosx smx , , 

Ms) = 3 — =- - 3 — , (2.13) 

X 6 X 1 X 

are the standard spherical Bessel functions. 

By substituting Eq. (U) into Eq. (|j|), and combining it with Eqs. (^9|)-(^T3D, one can 
derive the following expressions: 

4 

d^d {2)ab = J2 i A k cos(2Aa + A;A$ ) + B k sin(2Aa + kA%) , (2.14) 

fc=0 

4 

d^S b d (2)ac S c = Y, [Ck cos(2Aa + A;A$ ) + D k sin(2Aa + A;A$ )] , (2.15) 

k=0 
4 

c^S^S 6 ^ = ]T [S* cos(2A« + A;A$ ) + F fe sin(2Aa + kA$ )] , (2.16) 

fc=0 

where Aa and A$ are given by Eq. ( |2.5| ), and A k , B k , C k , D k , E k and F k (for k = 0, .., 4) 
are numerical coefficients, that are given in Appendix A. Inserting Eqs. ( [2.15|) -( p.l6|) in 
Eq. ( |2.7| ), the overlap reduction function becomes: 

4 

7(/) = E t p fe( x ) cos^(Aa ,A$) +Q fc (x) sin cp k {Aa , A$)] (2.17) 

fc=0 



where 



Pk{x) = pi(x) A k + p 2 (x) C k + p 3 (x) E k , 

Qk{x) = Pi(x) B k + p 2 {x) D k + p 3 {x) F k , (2.18) 
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FIG. 2. The overlap reduction function j(f) for a pair of identical laser interferometers 
with the LISA orbital configuration. The plot shows j(f) as a function of the frequency for four 
different centre-of-mass separation angles A<3?o, which correspond to different distances of the two 
detectors, cfr. Eq. ( fOll ): A<J» = 20° (solid line), A$ = 40° (dashed line), A<J» = 60° (dotted 
line), A$o = 90° (dotted-dashed line). Here we always assume that the initial orientation of the 
detectors is such that Aao = 0. The small panel at the top zooms the behavior of 7(/) in the mHz 
region. 

depend only on the detector separation and the radiation frequency, and 




(2.19) 



is a function of the relative orientation of the instruments. 
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FIG. 3. The overlap reduction function for a pair of identical laser interferometers with the 
LISA orbital configuration. The plot shows j(f) as a function of centre-of-mass separation angle 
A<I>o (in units 2ir), for selected frequencies: / = 10 -2 Hz (solid line), / = 10~ 3 Hz (dotted line), and 
/ = 10~ 4 Hz (dotted-dashed line). The constants «oi and ao2 are selected in such a way that, at 
the (arbitrary) reference time t = 0, Aao = 0. 



We would like to stress that our definition is such that 7(/) = 1 V/, for co-aligned and 
co-located interferometers with arms perpendicular to each others, that is the angle between 
m and n is ir/2. However, the LISA opening angle is 7r/3, and the detector response is 
reduced by the factor Vo/2: as a consequence, the maximum value that 'j(f) can attain is 
3/4. 

Notice also that, as pointed out by Cutler |30|, the read-outs from the three arms of 
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LISA can be combined in such a way to form the outputs, say oi and oji, of two co- located 
interferometers rotated by 7r/4 one with respect to the other, whose noise is uncorrelated at 
all frequencies. Unfortunately, the cross-correlation of oj with on is useless for searching for 
stochastic backgrounds, as the overlap reduction function is identically zero over the whole 
frequency range. 

Figs. and § show the behavior of 7(/) as a function of the frequency and the separation 
angle A$o, which is equivalent to the distance D, see Eq. ( |2.11p . Here we assume that the 
instruments are inserted into their orbits so that Acto = 0. In this case, at fixed frequency, 
the overlap reduction function depends only on A$ , as the relative orientation of the 



detectors is determined only by $(£), c f- Eqs. ( |2.3[ ) and (|2.4). This also implies that as the 
two instruments are moved apart - A$o increases - their orientation changes too. In order 
to achieve the maximum 7(/) for a given separation, one should therefore suitably tune 
the initial orientation of the detector arms. The plots clearly show that at low frequencies, 
/ ^ 5 x 10~ 4 Hz, 7(/) is fairly flat and close to its maximum value (which is set by the 
separation and orientation); in fact the radiation wavelength is A gw ~ 2(//lmHz) _1 AU, 
and for separations smaller than 1 AU, the degradation of SNR at very low-frequencies is 
less than a factor « 2. At high frequencies, say / £ Hz, placing two interferometers at 
a distance larger than 10 6 km, would severely degrade the SNR, by a factor of 10 or more. 



III. THE ASTROPHYSICALLY GENERATED STOCHASTIC BACKGROUND 

The so-called astrophysically generated GW stochastic background (GGB) is mainly due 
to the incoherent superposition of gravitational radiation emitted by short-period solar-mass 
binary systems. A variety of binary populations contribute to it, but the main contribution, 
in the mHz region, is due to close white-dwarf binaries (CWDB's). Present estimates suggest 
that it is above the LISA instrumental noise in the frequency region ~ 10" 4 - 3 x 10" 3 Hz, 
right at the heart of the observation window, see Fig. [TJ. However, sizeable effects are also 
given by other sources such as W UMa (Ursae Majoris) binary stars, unevolved binaries, 
cataclysmic binaries, neutron star-neutron star (NS-NS) binary systems, black hole-neutron 
star (BH-NS) binaries, and possibly BH-BH binary systems. 

The GGB is a guaranteed GW source in the low-frequency band; however, it is likely to 
overwhelm the PGB, degrading the sensitivity of the instruments in searching for a stochas- 
tic signal produced in the early Universe. A rigorous analysis of GGB's goes far beyond 
the purpose of this paper; here we review the main features and discuss the fundamental 



theoretical issues. We refer the reader to |pl|-p5|1, and references therein, for a thorough 
discussion of the astrophysical sources. 

Ultimately, the reason that the radiation generated by large populations of binary sys- 
tems is effectively a stochastic signal is simple: there are too many free parameters that one 
needs to fit the data in order to resolve all the binary systems that contribute to the signal. 
In fact, our galaxy contains ~ 10 7 CWDB's; they evolve, due to radiation reaction, over 
a time scale ^ 10 7 years. Therefore, during the typical observation time T ~ 1 yr, they 
are seen as highly monochromatic, and, in the band 10 -4 — 10~ 3 Hz, each frequency bin is 
"contaminated" by roughly 10 3 sources. The problem of resolving each individual binary 
system is actually made worse by the motion of the detector, because, in addition to the 
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frequency, one needs to solve also for the source position in the sky and the orientation of 
the orbital plane. 

In the following we discuss the estimates Q g (f) of the spectral energy density of generated 
backgrounds, and in particular their isotropic component, and the critical frequency f g up to 
which GGB's are indeed present. For / £ f g , the individual binary sources of astrophysical 
populations can be resolved, and their radiation subtracted from the data, opening up 
the band to search for a PGB. The former two properties are essential in addressing the 
possibility of detecting PGBs with space-based experiments. 

In this context, it is useful to divide the sources that contribute to the GGB into two 
categories: galactic sources and extra-galactic sources. Their key distinguishing feature of 
relevance here is the degree of isotropy of the GGB that they generate, for an observer 
on board a LISA-like detector. In fact, the extra-galactic contribution is expected to be 
isotropic to a rather high degree (the radiation being dominated by binary systems at cos- 
mological distances; for more details see f3TS||); it is, therefore, impossible to discriminate it 
from the PGB. On the contrary, the GGB produced by galactic sources is clearly highly 
anisotropic. In fact, galactic stars are spatially distributed, approximately, according to 
exp(— r/r ) exp[— (z/zo) 2 ], where r is the radial distance to the Galactic centre, r — 5 kpc, 
and z the height above the Galactic plane; Zq ~ 5 kpc for neutron star binaries and less 
than 300 pc for the other types. Due to the peripheral location of the Solar System, and 
the change of orientation of the LISA arms during the years-long observation time, galac- 



tic generated backgrounds appear strongly anisotropic fl37j . It is also conceivable that the 
isotropic portion of the galactic contribution does not exceed the total extra-galactic GGB. 

Several astrophysical uncertainties affect the estimates of the GGB that have been carried 
out so far. A careful analysis of the galactic contribution has been performed by Bender 
and collaborators , taking into account a wide range of binary populations. A good 

fit of the galactic GGB spectral density is: 

M Q-42.685 f- 1.9 l -5</<10- 3 - 15 , 

S g (f) = \ io- 60 - 325 /- 7 - 5 10- 3 - 15 < / < lO" 2 - 75 , (3.1) 

( 10-46-85 f-2.6 10 -2.75 < f ? 

which is related to £l g (f) by Eq. fll.flj). In Eq. (|3.1| ) the space density of CWDBs is assumed 
to be 10% of the theoretical value predicted by Webbink [38[|, and the radiation from helium 



cataclysmic variables, likely to contribute significantly in the frequency window ~ 1 mHz — 3 
mHz, is not taken into account. The estimate (|3.1|) can be definitely regarded as a solid 
lower limit, and is likely correct within to a factor ~ 3; see also |33| , |3^| . We would also like 
to stress that the dominant contribution from CWDBs switches off at ~ 10~ 2 Hz, where 
the binary systems coalesce; for / £ 10 -2 Hz, NS-NS binary systems are the sources that 
contribute most to the GGB. 

The extragalactic contribution to the GGB has been estimated in ||34j| , and is weaker than 



the galactic contribution by a factor 10-to-3 in the relevant frequency band, where the 
main uncertainty comes from the star formation rate at high redshifts [f3fJ,f|(J (see however 
the optimistic estimates in |35|| ). As we have mentioned at the beginning of the section, 
it is likely that the total (galactic + extra-galactic) isotropic component of the GGB does 
not exceed the total extra-galactic GGB; indeed, we will assume that the isotropic portion 
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of ftg(f), the only one that affects searches for the PGB, follows the frequency distribution 
predicted by Eq. (|3.1|) , but reduced by a factor e < 1. In the rest of the paper we will 
therefore assume 

4vr 2 

n 9 , ts = e—,fS(f). (3.2) 

Note that the determination of the value of e is a delicate matter. At frequencies above a 
few mHz, where the GWs from galactic sources can be subtracted, it is likely that e ^ 0.3, 
while below ~ 1 mHz, Ref. [57| suggests e ~ 0.7. However a more detailed analysis is needed 



in the future to provide a better estimate of this value. In the following we optimistically 
set e = 0.1, which can be regarded as a solid lower limit; the results that we will present in 
the next Section can be easily rescaled as a function of e. 

We consider now the critical frequency f g up to which radiation emitted by solar-mass 
binary systems (in the whole Universe) produce a GGB. For / ^ f g the observational window 
becomes "transparent" to the primordial GW background. Following the discussion at the 
beginning of the section, we estimate f g by using a very simple argument based on first 
principles: if the number of independent degrees of freedom of the data set - the number 
of data points, or, equivalently, the number of frequency bins - is smaller than the total 
number of independent parameters that describe the radiation, then the superposition of 
monochromatic GWs must be considered as a stochastic background. If the opposite is true, 
we have enough information to characterize, at least in principle, each individual source, and 
the signal is a deterministic one. For the sake of simplicity - although not exactly true, see 
the discussion at the beginning of the section and [|3(]] - we assume that each binary system is 
characterized by one parameter. The critical frequency f g is, therefore, formally determined 
by the condition that the (average) number of sources per frequency bin is less than one: 

^A b /<l; (3.3) 

here dN/df is the number of binary sources, emitting at frequency /, per unit frequency 
interval. Assuming that the merger rate is R, and in the relevant frequency range the 
binaries evolve only through radiation reaction, we have 



df 

where df /dt can be estimated using the Newtonian quadrupole formula: 

dt 5 

and M. = (mim2) 3//5 /(mi + 7712) 1 / 5 is the so-called chirp mass (mi and m2 are the masses 
of the two orbiting stars). Substituting Eqs. (|3.4j ), and (|3.5| ) into Eq. (|3.3|) , the frequency 
f g is easily determined: 

/ ^ \ 3 / n /R\ 3 / n 
/.= (!) ^ 8/11 A<- 5/11 (f) ■ (3.6) 
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We are now interested in determining an upper limit to f g , considering binary populations 
from the whole Universe. For / £ 10 mHz, the main contribution to the GGB is given by 



NS-NS binaries [ 31 1 . Their merger rate is uncertain; current estimates yield a galactic merger 



rate in the range [41-^3] : 



ins 



10" 6 -5 x Hr 4 yr _1 - (3.7) 



We can extrapolate this result to the entire Universe by simply multiplying the galactic rate 
by the total number of galaxies Nq'- 

By using this approach, we assume that i?NS does not vary with the redshift, which is 
probably not true. However, even if at high redshift i?NS is a factor 10 higher than in our 
galaxy, the very weak dependence of f g oc R 3 ^ 11 on the merger rate, see Eq. ( |3.6|) , ensures 
that this crude estimate is correct within a factor ~ 2. Assuming that the typical chirp mass 
of the binaries in the population is M. = 1.2 M , which corresponds to mi = m 2 = 1.4 M & , 
and the rate Q, Eq. (|3j) yields: 

This simple analysis leads therefore to the conclusion that the window / £ 0.1 Hz is likely 
free from stochastic backgrounds generated by astrophysical sources; radiation from NS-NS 
binaries is still present, but one can detect each individual source, estimate its parameters, 
and remove from the data stream the signals. In principle, the search for the PGB becomes 
limited only by the instrumental noise. 



IV. SENSITIVITY 

We can now proceed to discuss the sensitivity of LISA, and possible follow-up missions, 
to search for the PGB. The major disadvantage of the presently designed LISA mission is 
the lack of two instruments with uncorrelated noise; in fact, each pair of the three co- located 
LISA interferometers share one arm, and therefore common noise. One can, in principle, 
extract two data streams with uncorrelated noise at all frequencies |30| . unfortunately they 



are equivalent to the outputs of a pair of detectors, one rotated by 45° with respect to the 
other, and as we have stressed in Sec. |H], the response to a stochastic signal is then null. 
Nonetheless, LISA can play an important role in searching for stochastic backgrounds, as it 
can possibly achieve a sensitivity hl 00 Q ~ 10~ n — 10~ 10 by exploiting some intrinsic proper- 
ties of the signal and/or the unique features of the instrument to operate in a configuration 
where the GW signal is (almost) absent. 

During the observation time, LISA changes orientation with period 1 yr. A background 
which is anisotropic produces a periodic modulation in the auto-correlation function that 
can stand above the noise [|37|- Unfortunately, one can not exploit this feature to detect 
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a primordial signal that we expect to be intrinsically isotropic to a high degree; the only 
sizeable anisotropy that one can foresee is a dipolar one, produced by the motion of our 
local system (and therefore LISA) with respect to the cosmological rest frame with velocity 
^prop/c — 10~ 3 . Unfortunately, an interferometer has a quadrupolar antenna pattern, and 
the SNR produced by the quadrupole component is reduced by a factor (v prop /c) 2 ~ 1CT 6 . 
Nonetheless, the use of the signature induced by anisotropies could lead to the detection of 
galactic generated backgrounds, as suggested in |57] . 

Time delay interferometry with multiple readouts |Q provides a way of suppressing by 
several orders of magnitude the GW contribution from the LISA output at frequencies below 
a few mHz, providing a shield to GW radiation. By exploiting this feature, one can calibrate 
the noise-only response of the instrument, and search for an excess power in the data stream 
that can be assigned to a signal of cosmic origin. 

LISA can therefore carry out searches for stochastic backgrounds with a single instrument 
at a very interesting sensitivity level. LISA is likely to detect the galactic astrophysically- 
generated background, and could set important upper-limits on the primordial contribution, 
that can rule out existing models. It also provides invaluable information regarding the 
astrophysically generated backgrounds that might turn out to be crucial in designing future 
follow-up missions devoted to GW cosmology. 

However, in order to reach the sensitivity h1 00 Q ~ 10~ 16 that we have set as a goal, 
two separated interferometers are essential. It is important to understand whether this 
target sensitivity is within the reach of near future space technology, and to discuss possi- 
ble fundamental limitations that might prevent us from achieving a detection at the level 
h\ m VL ~ 10 -16 . The science goals are of such importance that some suggestions have been 
already put forward as to how to introduce modifications to the currently envisaged LISA 
configuration in order to accommodate a pair of independent instruments, possibly with the 
capability of shifting at-will the centre of the sensitivity window from ~ 1CT 3 Hz to ~ 0.1 

We start by considering a pair of identical LISA detectors; we therefore assume present, 
or near future, technology. Then, we discuss "second generation" LISA detectors, specifically 
aimed at PGB searches. 



A. The sensitivity of two LISA interferometers 

We consider what the present LISA technology allow us to achieve by cross-correlating 
the outputs of two identical interferometers located at a distance D, cf. Eq. (|2.11|) ; this is 
equivalent to determining the minimum value of ti\ m Q that one is able to detect. No unique 
answer can be given to this question, as it depends, of course, on the frequency dependence 
of the true signal fi(/). Lacking any solid theoretical prediction, we choose a "maximum 
ignorance" approach, and assume that is constant over the relevant frequency range. 
This hypothesis is not unreasonable, because the frequency band over which the SNR builds 
up is fairly small, due to S(f) and 7(/), cf. Eqs. (|2.6| ), ( |2.17| ), and Figs. |I]and By setting 
£l(f) = const, and S^(f) = S^'(f) = S n (f), and solving Eq. ( |1.20| ) for fi, we can estimate 
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the minimum detectable value of the energy density content of the GW background^: 



,2 (min) K V 507r 



inn - Tl/2 



7 2 (/) 2 



(4.1) 



the constant K is related to the false alarm probability and the detection rate associated 
with the measurement of a background with energy f2( min ). For a false alarm probability of 
5%, and a detection rate of 95%, we have K ~ 3.76 (the sum of the false alarm probability 
and the detection rate need not to be 1, and these two quantities are totally subjective) ]2"9"[ . 
We therefore compute Eq. ( |4.1| ), with noise spectral density and overlap reduction function 
given by Eqs. (|2.6|) and ( |2.17|) , respectively. The results presented here are computed, for 
sake of simplicity, for the case Aa = 0, and the integration is carried out over the frequency 
band 10" 4 - 10" 2 Hz. 

It is interesting to analyze first how different frequency regions of the whole sensitivity 
window contribute to the total SNR. Fig. (| shows the fraction of the total signal-to-noise 
ratio that is accumulated per unit logarithmic frequency interval; it indicates clearly that 
the key frequency band is ~ 8 x 10~ 4 Hz — 5 x 10 -3 Hz. 

Fig. |5| shows the sensitivity to a generic GW stochastic background that one could in 
principle achieve for a time of observation T = 10 7 sec. It is straightforward to rescale 
these results for a different observation time, as SNR oc T 1 / 2 and /i 2 00 f2^ mm ^ oc T -1 / 2 , see 
Eqs. ( |1 . 1 2| ) and ( |4.1| ). It is remarkable that LISA is capable of measurements in the range 



5 x 10~ 14 £ hf 00 fi( mm ) <, 10~ 12 , which is about three orders of magnitude better than can be 
achieved (although in a different frequency regime) with Earth-based interferometers oper- 
ating in the " advanced" configuration. The latter experimental set-up requires considerable 
technological and possibly conceptual developments, and is expected to be implemented not 
before the end of the decade. The one-order-of-magnitude range in /i 2 00 
effect of the overlap reduction function, where the lower limit is for co-located instruments, 
and the upper limit for detectors at a distance D = 2 AU with Aao = 0, cf. also Fig. [|[ 

However, such instrumental sensitivity does not correspond to a comparable sensitivity 
to PGBs; in fact, the peak of the GGB spectrum, of the order ~ 10~ n for the isotropic 
component, is where the instrument is most sensitive, cf. Fig. |l[ In order to quantify this 
effect, and to address the capabilities of space-based interferometers for cosmology, we need 
a short digression. 

Consider a generic two-component stochastic signal: 

n(/) = n a (/) + na(/), (4.2) 

and assume that ^i(/) and f^/) share exactly the same statistical properties (they are 
isotropic, stationary, Gaussian and unpolarized), and therefore can not be distinguished 
from each other. In order to search for Qi(f), one constructs the optimal filter 



1 Notice that it is appropriate to use the weak (with respect to the noise) signal approximation, 
Eq. ( 1.20| ), as we are aiming at the detection of the weakest possible background, which is clearly 
dominated by the noise. 
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FIG. 4. The contribution of different frequency regions to the total signal-to-noise. The 
plot shows the fraction of SNR 2 , Eq. ( |1.20| ), that is accumulated per unit logarithmic frequency 
interval, d(SNR 2 )/(f(ln/), as a function of the frequency /, for cross-correlations involving a pair of 
LISA instruments. The vertical axis is normalized so that SNR 2 = 1 for co-located and co-aligned 
interferometers, with an arm opening angle of 60°. The solid and dotted-dashed lines refer to 
D = and D = 1 AU, respectively. In both cases A«o = 0, and the noise spectral density is given 
by Eq. Q. 



and correlates it against the data of two instruments, following the scheme described in the 
Introduction, Sec. |TC| . The signal-to-noise ratio at the output of the filtering process is 
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FIG. 5. The minimum detectable value of the fractional energy density of a GW stochastic 
background. The plot shows h\ QQ Jl( mm ) - assumed to be flat, see Eq. Q4.1| ) - for a cross-correlation 
experiment involving two identical LISA detectors, as a function of the instrument separation D 
(in AU). The integration time is set to T = 10 7 sec. The instrument noise spectral density and the 



overlap reduction function are given by Eqs. (|2.6[) and (2.17), respectively; for the sake of simplicity, 
we set Aao = and the integration band corresponds to 10~ 4 Hz — 10~ 2 Hz. The dashed line refers 
to an experiment limited only by instrumental noise, where the false alarm probability and the 
detection rate are 5% and 95%, respectively; in this case h 2 100 n^ scales as T~ l l 2 . The solid 
line refers to an experiment limited by instrumental and confusion noise, and shows the minimum 
detectable value of the energy spectrum of a primordial background hi 00 Qp m ' . In this case, the 
spectrum of the astrophysical generated background is computed according to Eq. (|3^), and we 
have optimistically set e = 0.1. The limit on /i 2 00 fip mm ^ does not improve, in this case, by increasing 
the integration time and/or lowering the instrumental noise. We refer the reader to the text for 
further details. 



SNR = T 



IOtt 2 



f 3 R(f) 



Qi, Q 
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Q 1 ,Q 1 ) 1+ ) A A ( } . (4.4) 



The second term in brackets can be interpreted as the (undesired) residual correlation in 
the detection filter due to 0, 2 (f), which can not be eliminated. When it exceeds unity, the 
component can not be detected, regardless of the instrumental sensitivity; observations 
carried out with smaller noise disturbances, would simply increase both terms (Qi,Qi) 
and (Qi,Q2) by exactly the same amount, without improving the chances of detecting f^. 
The same holds for the integration time T: it has no effect at all on the capability of 
discriminating the two components. Therefore, the minimum detectable value of is 
set by the condition 



(Ql,Ql 



QxM ■ (4.5) 



Notice that the frequency dependence of the two components is very important: if f2i(/) and 
^2(/) follow a similar frequency behavior, the filter picks up more power from the "spurious" 
component Q^if)', if they are drastically different, even if fyi(f) dominates Q\(f), one could 
achieve a detection. 

This example describes exactly the issue that we are considering in this section, by 
simply identifying Q± with Q p , and 2 with Q gt i s . The unresolved radiation from binary 
systems provides therefore a fundamental sensitivity limit in searching for the primordial GW 



background. We have computed this limit from Eq. ( f4.5|) , in the case of an experiment carried 
out with a pair of identical LISA's, assuming Q p constant, and fl g ^ s (f) given by Eq. ( [3.21) . 
The results are summarized in Fig. ||; the key conclusion is that two LISA detectors will be 
able to detect a PGB (with constant energy spectrum) only if hf 00 Q p £ 5 x 1CT 13 . The big 
loss of SNR with respect to the case where the experiment is limited only by instrumental 
noise is due to the fact that the GGB is very strong in the mHz band, and fl p (f) and Q g (f) 
have a similar decreasing frequency behavior in the frequency window where most of the 
SNR is accumulated: the residual correlation at the filter output produced by the GGB is 
large. 

Given the results reported in Fig. ||, it is straightforward to conclude that the GGB is 
a guaranteed, strong GW signal for space-based detectors. In fact, if one constructs a filter 
matched to a GGB given by Eq. ( |3.2|) , and performs cross-correlations between two identical 
LISA instruments which are characterized by the noise curve (|2.6|) , one can detect such a 
signal with SNR g ~ 100, for a time of observation T = 10 7 sec; see Fig. |!J two to three 
days of integration time are sufficient to reach SNR g ~ 10. Two LISA-like detectors would 
therefore be extremely powerful telescopes to launch deep surveys of populations of binary 
systems in our Universe with periods between a few hours and a few hundred seconds. 

The bottom line of this analysis is therefore clear: the fundamental limiting factor in 
searching for a primordial GW background in the 10 -5 Hz — 10 -2 Hz frequency window is the 
stochastic radiation from unresolved binary systems. Searches for PGBs with h\ m Q p £ 5 x 
10~ 13 thus call for a change in the observational window. We briefly discuss this issue in the 
next Section. 
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FIG. 6. The signal-to-noise ratio at which the astrophysically generated background can be 
detected. The plot shows the SNR that can be achieved in 4 months of integration time, correlating 
the data of a pair of LISA detectors, as a function of the distance D (in AU). The Wiener filter is 
matched to the signal described by Eq. (|3.2|). 



B. Towards testing slow-roll inflation 



"-100 "p 



The mHz frequency window is unsuitable to reach the ambitious sensitivity level 
10~ 16 predicted by slow roll inflation. One needs to design an experiment 
with optimal sensitivity in a band free from generated backgrounds. Given our present as- 
trophysical understanding, the most promising region seems to be ~ 0.1 Hz — 1 Hz, which 
would be optimally accessible through space-borne interferometers with arms shorter than 
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the LISA ones by a factor ~ 100. In fact, the entire frequency window from ~ 10~ 7 Hz to 
~ 0.1 Hz is contaminated by stochastic signals of astrophysical origin with /i 2 00 Q ^> 10~ 16 . 
A background generated by massive black hole binary systems with fractional energy density 
/z 2 00 Q ~ 10~ 15 — 10 -14 is present in the /zHz range f46fl ; however, our ignorance concerning 
the formation rate of massive black holes and the merger rate of massive black hole binary 
systems in the range ~ 10 5 M & — 10 9 M & prevents us from giving a more solid estimate of the 
background generated by such objects. For 10~ 5 Hz £ f £ 10~ 2 Hz, unevolved binaries and 
WD- WD binary systems completely swamp the observational window; see Sec. [HT[ Above 
~ 10 niHz the only residual sizeable contribution comes from NS-NS binary systems; as we 
have discussed in Sec. (HI], above ~ 0.1 Hz, the number of sources per frequency bin becomes 
less than one, and the sky is "transparent" to a primordial signal. For rather long integration 
times ~ 3 yrs, the rms instrumental noise level that is required to test the prediction from 
slow-roll inflationary models is of the order ~ 10 -24 
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Operating at considerably higher frequency than LISA, the two detectors would have to be 
closely located, D ^ 10 11 cm, in order to have optimal overlap reduction function; however, 
this also increases potentially correlated noise sources. Clearly the technological challenge to 
achieve the mentioned sensitivity is considerable; the main noise sources that would degrade 
the performance of such a detector are the shot noise, beam pointing fluctuations, and 
accuracy of the phase measurement technique. This imposes stringent requirements on the 
power and frequency of the laser, as well as on the dimensions of the "optics" and on other 
components of the instrument. We would also like to stress that the value of h vms that we 
have quoted in ( [4.6|) is the effective noise fluctuation in the data stream, after the spectral 
lines from individual NS-NS binaries, that are still copiously present in the observation 
band, have been removed. How this can be effectively done and what instrument sensitivity 
is required is still an open question that requires a careful analysis ||47 |. The main result of 
this crude analysis is that one could indeed reach the target sensitivity, and the fundamental 
limitations that make the mHz band unsuitable are removed. However, several questions 
remain open for future analysis: the discussion of these issues goes far beyond the purpose 
of the present paper. 



V. CONCLUSIONS 

Gravitational wave experiments in the low-frequency window, together with high- 
frequency ground-based interferometers, are expected to improve our picture of the very 
early Universe, and the understanding of the behavior of fundamental fields at high energy, 
by detecting, or setting stringent upper limits on the primordial background of gravitational 
radiation. 

In this paper we have analyzed the sensitivity of space-borne laser interferometers of 
the LISA class, and possible succeeding missions. In order to set a reference frame for this 
discussion, we have regarded the detection of a GW background produced during the early 
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Universe of energy h\ Q0 ft p ~ 10~ 16 , consistent with the prediction of standard slow-roll infla- 
tion, as the goal of GW cosmology. We have assumed the operation of two space-detectors, 
in order to achieve the best sensitivity and detection confidence, and we have shown that the 
technology available for LISA already ensures the detection of a GW background as weak 
as h1 00 Q ~5x ICC 14 . However, the strong stochastic signal in the mHz band due to short- 
period solar-mass binary systems that can not be resolved as individual sources prevents us 
from detecting a primordial background weaker than h\ m Q p « 5 x ICC 13 . Astrophysically 
generated stochastic backgrounds therefore set a fundamental limit in the mHz band that 
prevents us from achieving a sensitivity that goes beyond what is already guaranteed by the 
LISA technology. They also represent a guaranteed strong signal detectable at high signal- 
to-noise ratio, which enables the study of the distribution and merger rate of populations of 
binary compact objects in the Universe. 

Dedicated missions with optimal sensitivity in the window 0.1 Hz - 1 Hz appear, at 
present, the only viable option in the search for very weak primordial backgrounds, and 
we have briefly discussed the technological challenges involved in probing slow-roll inflation. 
Our order-of-magnitude analysis strengthens the hope that a sensitivity level h\ Q0 Q p ~ ICC 16 
might be within the capability of future dedicated low-frequency detectors. 

Our analysis clearly indicates the key issues that deserve further investigation: a solid 
estimate of galactic and extra-galactic GW backgrounds produced by astrophysical sources, 
the investigation of the statistical issues that can lead to the discrimination of the PGB 
from the GGB, and a more rigorous analysis of the technical and conceptual problems for 
low-frequency experiments dedicated to GW cosmology. On the observational side, the 
presently designed single-instrument LISA mission is a fundamental step for the planning of 
more ambitious, multi-detector experiments: we will be able to measure directly the degree 
of anisotropy of the generated background, shedding light on the fundamental limiting factor 
of mHz experiments. In fact, while the present paper deals only with the detection of an 
isotropic stochastic signal, the remarkable sensitivity of LISA offers the chance of going far 
beyond: a detailed study of the anisotropy and angular dependence of stochastic signals, 
both of astrophysical and primordial origin. Such an investigation is currently in progress, 
and will be reported in a separate publication p8 . 
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APPENDIX A: 



We give here the values of the numerical coefficients that enter expressions ( 2.15Q , ( 2.16|) 
and ( |2.16| ) of the overlap reduction function derived in Sec. ( [11 B| ): 



27 



513 135 243 9 33 

4096 ' 1 ~~ 1024 ' 2 ~~ 2048 ' 3 ~ ~ 1024 ' 4 ~~ 4096 ' 

27^3 _ 27^ 81^ _ 27^ 27^ 

~ 4096 ' 1 ~~ 1024 ' 2 ~ ~~ 2048 ' 3 ~~ 1024 ' 4 ~ ~ 4096 ' 

513 C7= iIL C = — C =— C 33 - 
8192 ' 1 2048 ' 2 4096 ' 3 2048 ' 4 8192 ' 
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